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Abstract 


This  investigation  compared  thin-airfoil  theory  predict¬ 
ions  of  change  in  zero- lift  angle  of  a. tack  divided  by  flap 
deflection  versus  flap-hinge  location  fE)  and  change 

in  moment  coefficient  about  the  aerodynamic  center  divided  by 
flap  deflection  versus  E  to  wind  tunnel  studies  of 

the  same  effect  for  a  cambered  airfoil.  This  comparison  for 
cambered  airfoils  has  not  been  previously  reported.  The  NACA 
4412  airfoil  was  chosen  for  this  study.  Results  of  this  study 
show  that  thin-airfoil  theory  does  predict  the  trend  of  the 
data.  The  Aulo/h  versus  E  prediction  favorably  matched  the 

f 

data.  The  AC^pac/n  versus  E  data  differed  from  the  predictions 
in  a  similar  way  to  that  reported  for  uncambered  airfoils? 
however,  the  data  in  this  investigation  showed  less  favorable 
agreement  to  thin-airfoil  theory  than  the  previously  reported 
results.  This  difference  might  be  explained  in  part  by  the 
experimental  setup  used  in  this  study. 
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WIND  TUNNEL  INVESTIGATION 
OF  VARYING  HINGED  FLAPS 


I  Introduction 


Obi'CCtive 

objective  of  this  thesis  was  to  compare  thin-airfoil 
theoretical  predictions  of  the  change^in  zero- lift  angle  of 

attack  divided  by  flap  deflection  (AaLQ/n)"  versus  flap-hinge 

f  ' 

location'  (E)i‘and  the  change  in  moment  coefficient  about  the 
aerodynamic  center^ivided  by  flap  deflection  (AC^^c/d) 
versus  E  to  wind-tunn)^  studies  of  the  same  effects  for  the 
moderately  cambered  NACA  4412  airfoil. 


Background 

In  thin-airfoil  theory,  the  airfoil  is  replaced  with  its 
mean  camber  line.  Point  vortices  with  a  differential  strength 
of  ydx  are  then  distributed  continuously  along  the  mean  camber 
line  to  induce  velocities  which,  when  added  vectorially  to  the 
free-stream  velocity,  produce  velocities  tangential  to  the  mean 
camber  line.  Because  the  maximum  camber  is  usually  small  in 
comparison  to  the  chord,  the  boxindary  conditions  for  the  flow 
can  be  satisfied  on  the  chord  Instead  of  at  the  mean  camber 
line.  This  produces  an  induced  velocity  (V^n)  x©  such  that 


(1) 


1 


at  each  point  on  the  chord  line.  With  the  assumption  of  small 
angles,  the  free- stream  velocity's  normal  component  to  the 
mean  camber  line  (V^^)  is  such  that 

V«.n  =  V„[a-^]  (2) 

The  induced  velocity  must  also  be  equal  to  the  negative  of 
at  Xq  if  there  is  to  be  no  flow  through  the  airfoil,  i.e. 

.  So,  from  equations  1  and  2 

-  ''.[i'-o)  -  a]  <3, 

By  specifying  a  y(x) -distribution,  the  shape  of  the  mean 
camber  line  and  a  to  produce  such  a  distribution  can  be  found 
or  vice  versa.  The  only  restriction  on  the  y(x) "distribution 
is  the  Kutta  condition,  y(c)  =  0  ,  which  must  be  applied  to 

determine  the  circulation  and  keep  Vj^^(c)  finite  so  that  the 
flow  can  leave  smoothly  and  tangentially  from  the  trailing  edge 
as  in  real  viscous  flow. 

By  using  the  transformation  x  =  ~  cos  0)  ,  the  relation¬ 

ship  [cos  (n-l)0  -  cos  (n+1) 0] /2  =  sinn0sin0  ,  the  definite 

-  "[^]  .  T  into 

a  fourier  series  of  the  form 

^  f  A„slnne]  (A) 

^  n=l  J 

the  induced  velocity  or  left  hand  side  of  equation  3  now 
becomes 


Vln(«)  -  V.[-A„  +  f^A„cosn6j 


(5) 


Equations  3  and  5  imply  that 

dz 


^  “  (o  -  Aq)  +  J^AnCOS  ne 
n=l 


(6) 


Equation  6  has  the  form  of  the  cosine  series  expansion  of 
For  a  given  mean  camber  line  ^  is  a  known  function  of  0 
and  therefore  the  values  of  Aq  and  Aj^  may  be  written  directly 


as 


and 


dz 

Hi 


d0 


(7) 


V-l/o 


dz 

Hx 


cos  n0d0 


(8) 


Since  the  section  lift(8.)  equals  the  density  times  the 
free- stream  velocity  times  the  circulation  (pV„r)  where  r 
equals  the  integral  of  differential  y  distribution  over  the 
chord  U'  ydj^  ,  the  aerodynamic  coefficient,  Cj^,  defined  as 


£/-jpV^c,  can  now  be  rewritten  as 


2  r"" 

^JQ 


(9) 


By  using  equation  4  for  y,  and  carrying  out  the  integral  in 
equation  9,  Cg^  can  be  found  in  terms  of  A^  and  A^  as 


Cg  -  2irAj,  +  wAj 


(10) 


Replacing  Aq  and  A^  with  their  integral  forms  from  equations 
7  and  8 

Cg  -  2i»  [-i/o  ||(cos  0  -  Ddej  (11) 
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The  angle  of  attack  which  makes  the  lift  or  the  lift  coef¬ 
ficient  in  equation  11  equal  to  zero  is  called  the  angle  of 
zero  lift,  «lq*  Applying  the  definition  of  a^Q  to  equation 
11  gives 

“Lo -irTo  at  ”2) 

A  similar  development  for  will  give 

TT 

Cmac  =  iX  ll 

(Refs  1;AA-A7  and  2:116-132) 

Equations  12  and  13  are  used  in  Section  II  to  predict 
influences  from  varying  sizes  of  simple  hinged  flaps. 


A 


c. 


II  Numerical  Solution  for  Thin-Airfoil 


Cambered  Airfoils 

As  we  saw  in  Section  I,  thin-airfoil  theory  produces 
integral  equations  for  and  ^lq-  By  assuming  that  the 

camber  line  is  made  up  of  straight  line  segments,  these 
equations  can  be  rewritten  as  follows  (Ref  3:34) 

“LO  "  ^  e  - 1)  de  =  e  -  1)  de 

(14) 


and 


mac 


"  ifo 


^  (cos  20  -  cos 


0)  de  -  ^ 

i=l>-“*J  i-'e 


®i+l 


(cos  26  -  cos  6)  d0 


(15) 


or,  after  evaluating  the  straight  line  segment  integrals 


“LO  -  ^  H+l  -  0i+l>  -  (=1"  "i  -  8l)]] 


(16) 


and 


C 
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(17) 


These  equations  can  be  numerically  integrated  through  the 

dz 

summation  sign  by  using  a  digital  computer  as  long  as 

and  determined  for  each  straight  line  segment. 

The  correctness  and  accuracy  of  the  above  equations  can 
be  checked  by  evaluating  and  Cmac  analytically  and  then 
numerically  with  the  computer  for  the  same  airfoil.  Since 
this  Investigation  used  the  NACA  4412  airfoil  which  has  its 
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camber  line  in  functional  form,  both  methods  can  be  accomplish¬ 
ed  for  the  zero  flap  deflection  case.  The  analytic  solutions 
as  given  in  Appendix  A  yield  =  -4.1544826“  and  = 

-.1062391.  The  numerical  computer  solution  using  400  straight 
line  segments  yields  a^Q  =  -4.1540554“  and  =  -.1062276. 

While  the  computer  solution  is  a  crude,  brute  force  method,  it 
gives  extremely  accurate  results.  Appendix  B  contains  the 
program  from  which  and  for  the  4412  were  computed. 

Flapped  Airfoils 

The  influence  of  small  flap  deflections  on  a.  and  C 

LO  mac 

can  also  be  found  through  computer  numerical  integration. 

Figure  la  shows  the  coordinate  system  arrangement  for  an 
arbitrarily  cambered  thin  airfoil  at  an  arbitrary  angle  of 
attack  under  which  the  a|^Q  and  numerical  equations  hold. 

By  choosing  any  point  along  the  cambered  line  as  a  hinge  point 
for  a  flap  deflected  by  an  amount  n,  the  thin  airfoil  in  Fig. 
la  changes  to  that  in  Fig.  lb.  olq  may  be  computed 

for  the  flap  deflected  airfoil  by  defining  a  new  coordinate 
system  (z'  and  x')  as  shown  in  Fig.  Ic.  Before  the  a^o  and 
C^narS  of  the  flapped  and  unflapped  airfoils  can  be  compared 
to  arrive  at  Aolo/’I  and  for  particular  flap-chord 

ratios,  corrections,  if  they  exist,  for  using  the  z'  and  x' 
coordinates  as  the  computational  axes  must  be  applied  to  the 
new  and  Because  is  an  angle  dependent  upon  the 

coordinate  system  for  its  value,  the  angle  (a'  -  o)  between 
the  coordinate  systems  in  Fig.  Ic  must  be  subtracted  from  the 
new  or  flapped  airfoil  o^q.  No  correction  is  necessary  for 


6 


z 


Fig.  lb  Flap  arbitrarily  chosen  and  deflected. 


Fig.  Ic  New  coordinate  system  (x'.z')  to  find  new 


a-tn  and 
LO  mac 


7 


the  new  C^ac  since  it  is  related  to  the  airfoil  and  not  the 
coordinate  system.  By  doing  this  for  many  points  along  the 
camber  line  of  an  NACA  4412  airfoil  one  can  obtain  the  curves 
in  Figs.  2  and  3  depicting  the  influence  of  small  flap  de¬ 
flections  on  and  Reference  Appendix  B  for  the  pro¬ 

gram  used  to  find  the  data  points  for  Figs.  2  and  3. 

Figures  2  and  3  show  flap  influences  for  a  4Z  cambered  air¬ 
foil.  If  the  camber  is  different,  then  slightly  different 

curves  will  result  for  AC  versus  E  with  the  differences 

mac 

for  Auj^q/ti  versus  E  due  to  camber  being  negligible.  This  effect 
of  camber  on  the  AC^nac/’^  can  be  shown  by  comparing  4%,  IX,  and 
uncambered  airfoils  as  in  Fig.  4.  In  general,  even  the  largest 
difference  among  the  cambered  airfoils  and  uncambered  airfoil 
is  small.  The  above  results  can  be  checked  by  comparing  them 
to  predictions  which  have  not  been  derived  in  a  nijmerical  man¬ 
ner.  Chow's  (Ref  2:133-137)  analytic  derivation  of  Aulq/h 
versus  E  and  AC^^c/^  versus  E  for  uncambered  airfoils  results 
in  exactly  the  same  curves  derived  through  numerical  integra¬ 
tion  for  uncambered  airfoils,  as  they  should  be  in  order  for 
the  numerical  integration  approach  to  be  valid.  Chow  goes  on 
to  state  that  because  all  angles  are  small,  it  is  sufficient 
to  find  the  properties  of  an  uncambered  airfoil  at  zero  angle 
of  attack  with  flap  deflected  and  add  them  directly  to  the  prop¬ 
erties  of  the  cambered  airfoil  at  any  angle  of  attack.  That 
statement  agrees  with  the  small  differences  noted  between  the 
xmcambered  and  cambered  airfoils  for  ACmg^/n  versus  E  and  also 
serves  to  validate  the  numerical  approach  as  explained  here. 


Fig.  2  Effect  of  flap-chord  ratix>, 


Pig.  3  Effect  of  flap-chord  ratio  on  AC 


Effect  of  camber  on  versus 

mac 


Ill  Wind  Tunnel  Investigation 


Models 

The  models  for  wind  tunnel  testing  were  constructed  by 
the  AFIT  Model  Fabrication  Division.  The  basic  components 
of  a  test  model  are  shown  in  Fig.  5.  Components  common  to 
all  models  were  the  flap  deflection  mechanism,  the  side 
plates,  and  the  rear  lift  bar.  The  airfoil  component  con¬ 
sisted  of  two  sections  that  varied  because  of  changes  in  the 
flap-hinge  location.  The  change  in  flap -hinge  location 
thereby  produced  seven  different  models  with  flap-hinge 
locations  at  20Z,  30Z,  AOZ,  50Z,  60Z,  70Z,  and  80Z  of  chord, 
or  expressed  in  flap-chord  ratio,  .2,  .3,  .4,  .5,  .6,  .7,  and 
.8,  respectively. 


rear  lift  bar 


Fig.  5  Basic  components  of  a  fully  assembled  model 
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The  flap  deflection  mechanism  consisted  of  two  identical 
brass  plate.  see  Fig.  6.  Holes  B  through  G  served  to  attach 
the  brass  plates  to  cnly  the  flap  section  of  the  airfoil  be¬ 
ing  tested.  Hole  A  of  the  brass  plates  allowed  headed  screws 
to  be  freely  inserted,  through  both  the  side  and  brass  plates, 
into  treaded  holes  in  the  airfoil  flap  section.  This  estab¬ 
lished  an  axis  of  rotation  about  which  the  flap  section  was 
moved.  Holes  1  through  7  of  the  brass  plates  made  it  possible 
to  set  the  flap  at  prescribed  angles  of  deflection  when  used 
in  conjunction  with  corresponding  holes  in  the  side  plates. 
Appendix  D  gives  a  detailed  explanation  of  how  the  brass 
plates  were  used  to  establish  a  desired  flap  setting. 


Fig.  Brass  flap  deflection  plate  with  flap 

attachment  holes  (B-G) ,  flap  deflection 
holes  (1-7),  and  flap-hinge  hole  (A). 

The  two  side  plates  brought  all  the  components  together 
to  form  a  fully  assembled  model  ready  for  testing.  The  side 
plates  were  constructed  of  3/8-inch  clear  plexiglass  with 
beveled  edges  and  rounded  comers.  These  plates  served  three 
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additional  functions.  First,  they  served  as  end  plates  to 
simulate  two-dimensional  flow.  Secondly,  they  captured  the 
airflow  far  enough  upstream  of  the  airfoil's  influence  so 
that  the  tunnel's  free  stream  value  for  the  dynamic  pressure 
could  be  used  in  the  two  dimensional  analysis.  Finally,  the 
light  weight  and  transparency  of  the  plexiglass  side  plates 
facilitated  the  change  from  one  model  to  another  and  allowed 
viewing  of  the  model  during  testing. 

The  rear  lift  bar  was  constructed  from  tubular  aluminum 
with  the  cross  section  shown  in  Fig.  7.  This  bar  was 
connected  to  the  rear  lift  wire  of  the  wind  tunnel.  The  bar 
design  for  rear  lift  wire  attachment  eliminated  possible 
interference  problems  between  the  deflecting  flap  and  the 
rear  wire  attachment  mechanism.  It  did  minimize  the  number 
of  model  measurements  needed  to  set  angles  of  attack  and/or 
flap  deflections.  It  also  kept  tunnel  calibrations  to  a 
minimum  since  only  one  in-tunnel  calibration  of  a  model  was 
required  to  be  able  to  set  any  angle  of  attack  and  flap 
deflection  combination.  Additionally,  it  eliminated  any 
construction  or  structural  problems  associated  with  having 
the  rear  lift  attachment  mechanism  as  part  of  the  airfoil. 
While  the  bar  eliminated  many  possible  problems  and  minimized 
many  operational  problems,  it  was  not  wholly  satisfactory. 

The  shape  of  the  bar  as  shown  in  Fig.  7  led  to  data  errors 
which  had  to  be  accounted  for  as  discussed  in  Section  IV. 

The  seven  NACA  4412  airfoils  used  for  the  investigation 
were  made  of  solid  aluminxim  based  epoxy  resin.  Each  airfoil 
consisted  of  two  matched  but  separate  sections:  a  front 


Fig.  7  Cross  section  of  lift  bar  drawn  to  actual 
size. 

section  that  remained  fixed,  and  a  flap  section  which  rotated 
relative  to  the  side  plates  and  the  front  section.  Figure  8 
shows  the  seven  airfoils  and  the  difference  in  size  of  the 
flap  section  among  airfoils.  Each  front  section  had  two 
attachment  fittings  in  the  leading  edge  for  the  drag  and  front 
lift  wires.  The  matched  airfoil  sections,  which  formed  4A12s 
when  put  together  at  zero  flap  deflection,  were  made  from 
different  one  piece  4412s  because  of  the  chord  distance  and 
front  section  surface  lost  when  a  flap  section  was  cut  and 
shaped.  Thus,  a  flap  section  of  E  =  .6  was  matched  with  an 
appropriately  shaped  front  section  that  resulted  from  produc¬ 
ing  a  flap  section  with  at  least  a  flap-chord  ratio  of  .1 
less  than  .6.  The  one  piece  4412s  came  from  a  mold  made 
from  a  wooden  4412  airfoil  of  two  foot  span  and  nine  inch 
chord.  The  mold  insured  that  each  model  was  the  same  for  all 
practical  purposes.  The  similarity  of  models  was  checked 
through  model  measurement  prior  to  testing.  It  was  confirmed 
by  the  excellent  agreement  between  models  for  a^o 
slope  of  the  versus  a  curve. 

See  Appendix  C  for  detailed  model  geometry  and  meastnre- 


ments . 


Wind  Tunnel 

The  AFIT  Five-Foot  Wind  Tunnel  was  used  for  the  wind  tun¬ 
nel  investigation.  It  is  an  open  return,  closed  test  section 
wind  tunnel  with  a  5  ft  circular  cross  section  as  shown  in 
Fig.  9.  The  two  counter  rotating  12  ft  fans  of  the  wind  tun¬ 
nel  are  powered  by  four  AOO-horsepower  direct  current  motors 
easily  capable  of  producing  the  tunnel's  top  empty  speed  of 
350  mph.  The  speed  limit  prevents  further  cracking  of  the 
fan  blades  at  high  power  settings.  The  aerodynamic  forces  on 
the  models  were  determined  by  using  the  wind  tunnel's  wire 
balance  measuring  system  consisting  of  three  springless  scales 
connected  to  the  model  through  separate  wire  systems.  The 
wind  tunnel's  wire  balance  system  was  checked  and  aligned 
prior  to  beginning  the  runs  for  each  model. 

Installation 

Prior  to  installation  in  the  wind  tunnel,  the  holes  in  the 
side  plates  were  covered  with  tape  on  the  sides  facing  the 
t\mnel  wall  to  prevent  3-D  effects  and  the  gaps  between  the 
front  airfoil  section  and  the  flap  section  were  sealed  with 
modeling  clay.  Each  model  was  then  carefully  measured  on  a 
surface  plate  using  an  inclinometer  and  a  height  gage.  The 
measurements  were  used  to  determine  the  angle  of  balance  for 
the  wire  system.  The  angle  of  balance  then  determined  the 
distances  necessary  for  setting  angles  of  attack  for  the 
model.  Appendix  C  contains  a  more  detailed  explanation  of  the 
model  geometry  and  calculations  Involved. 

The  model  was  then  Installed  In  the  test  section.  Instal- 
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AFIT  Five-Foot  Wind  Tunnel  setup. 


lation  was  accomplished  by  connecting  the  model's  attachment 
fittings  to  plates  permanently  connected  to  the  appropriate 
tunnel  wire  system.  These  connections  allowed  the  model  to 
rotate  freely  at  each  connection  and  thereby  accurately  trans¬ 
mit  the  aerodynamic  forces  produced  on  the  model  to  the  wire 
system  and  then  to  the  scales  for  a  direct  measurement.  The 
two  trunnions  on  the  leading  edge  were  connected  to  the  drag 
and  front  lift  wires  while  the  attachment  fitting  on  the  rear 
lift  bar  was  connected  to  the  rear  lift  wire.  Figures  10a 
and  10b  show  a  model  installed  in  the  wind  tunnel. 

After  Installing  the  model,  the  wind  tunnel  instruments 
used  to  set  the  model  at  a  particular  angle  of  attack  were 
carefully  calibrated  prior  to  taking  data.  The  calibration 
used  an  inclinometer  to  measure  accurately  the  model's  angle 
of  attack  with  respect  to  the  tunnel  axis. 

Investigation  Procedure 

After  calibrating  an  Installed  model,  it  was  ready  for 
testing.  The  investigation  objective  required  that  enough 
data  be  gathered  to  plot  both  the  change  in  zero-lift  angle 
of  attack  divided  by  flap  deflection  versus  flap- 

hinge  location  (E)  and  the  change  in  the  moment  coefficient 
about  the  aerodynamic  center  divided  by  flap  deflection 
(ACmar/n)  versus  E.  Data  for  Aolq/^  versus  E  was  obtained 
by  first  setting  the  flap  at  0* ,  4* ,  or  8® ,  and  the  taking 
data  for  angles  of  attack  of  -6®,  -4®,  -2®,  0®,  2®,  4®,  6®, 

8®,  and  10®.  Data  for  AC^Qac/n  versus  E  was  obtained  by  hav¬ 
ing  the  model  set  at  0*  angle  of  attack  with  the  flap 


t. 


deflected  at  settings  from  0“  to  18®.  Overlap  in  gathering 
the  data  occurred  as  shown  in  Table  I  since  data  for  Aolq/ti 
versus  E  at  0°  angle  of  attack  was  also  used  as  data  for 
ACjuac/^  versus  E. 

A  typical  sequence  of  events  used  to  gather  data  for  a 
model  is  given  in  Table  I.  A  typical  change  from  one  flap 
setting  to  another,  an  event  in  Table  I,  obtained  data  for 
ACmac/’^  went  as  follows:  The  model's  flap 

setting  was  changed  by  one  degree  which  also  required  re¬ 
sealing  the  gaps  between  the  front  airfoil  section  and  flap 
section;  and  taping  over  any  uncovered  holes  in  the  side 
plates.  Next,  the  model  and  tunnel  were  inspected  to  insure 
no  loose  objects  remained  in  the  tunnel.  The  tunnel  doors 
were  then  secured,  after  which  the  tunnel  operator  aligned 
the  model  and  took  three  before-static  readings.  He  also 
ensured  a  current  barometer  and  temperature  reading  were 
taken.  The  tunnel  air  was  turned  on  and  brought  up  to  the 
testing  dynamic  pressure  of  three  inches  of  water.  The 
operator  adjusted  the  model  as  the  airflow  increased  to  keep 
it  properly  aligned.  With  the  model  properly  aligned  and  the 
dynamic  pressure  at  three  inches,  the  operator  took  five  wind- 
on  readings.  The  tunnel  air  was  then  turned  off  and  the  flow 
allowed  to  go  to  zero.  The  tunnel  operator  again  ensured  the 
model  was  aligned  properly  and  took  three  after-static  read¬ 
ings.  The  process  for  an  event  as  listed  in  Table  I  was  then 
complete  and  the  tunnel  could  again  be  entered  to  change  the 
model's  flap  setting.  Data  for  Aax,o/^  also  required  changes 
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TABLE  I 


in  angle  of  attack  for  the  flap  settings  of  0®,  4®,  and  8®. 
The  changes  In  angle  of  attack  were  accomplished  from  outside 
the  test  section  and  allowed  the  tunnel  air  to  remain  on  from 
one  angle  of  attack  to  another.  Five  wind-on  readings  were 
taken  for  each  angle  of  attack  while  three  static  readings 
were  taken  for  each  angle  of  attack  during  the  before-statics 
and  after-statics  portions  of  the  event. 

Tuft  studies  were  conducted  on  the  .2,  .4,  .5,  .6,  and  .8 
E  models  after  obtaining  the  force  data  for  that  individual 
model.  These  visual  studies  indicated  the  extent  of  separa¬ 
tion  on  the  airfoil's  upper  surface  and  were  used  later  to 
eliminate  data  for  those  particular  angles  of  attack  and 
flap  settings  which  showed  separation.  Appendix  E  gives  the 
results  of  the  tuft  studies. 

Data 

The  data  readings  for  any  one  setting  of  angle  of  attack 
and  flap  deflection  usually  consisted  of  three  before- static 
readings,  five  wind-on  readings,  and  three  after-static 
readings.  The  multiple  readings  helped  to  identify  bad  data 
and  reduce  scatter  through  averaging.  Additionally,  repeat 
readings  were  taken  to  verify  the  consistency  and  repeatabil¬ 
ity  of  the  data.  The  static  readings  determined  the  loads 
on  the  model  when  the  airflow  was  zero  while  the  wlnd-on 
readings  determined  the  loads  with  air  flowing  over  the  model. 

Each  time  a  reading  was  taken,  the  respective  loads  on 
all  three  scales  (drag,  rear  lift,  and  front  lift)  were 
printed  on  separate  paper  tapes.  The  readings  on  each  tape 


were  transferred  to  respective  tables  with  line  entries  con¬ 
taining  static  and  wind-on  readings  gathered  for  a  particular 
angle  of  attack  and  flap  setting  combination.  An  average 
static  value  for  any  combination  was  obtained  by  averaging 
all  the  static  readings  for  that  combination.  An  average 
wind-on  reading  consisted  of  only  averaging  the  five  wind-on 
readings  for  a  particular  event.  In  other  words,  any  set  of 
five  wind-on  readings  were  averaged  and  treated  as  a  separate 
result  from  any  other  five  wind-on  readings.  The  aerodynamic 
forces  on  the  model  for  any  particular  combination  of  angle 
of  attack  and  flap  setting  were  found  by  subtracting  the 
wind-on  averages  from  the  static  averages  to  give  drag  forces, 
rear  lift  forces,  and  front  lift  forces.  Appendix  F  contains 
the  raw  data  (data  attainable  directly  from  the  scale  tapes) 
obtained  during  the  investigation  for  each  of  the  seven 
models . 
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IV  Results 


«• 


Data  Reduction 

The  data  in  Appendix  F  was  first  reduced  to  obtain  lift 
coefficient  versus  angle  of  attack  graphs  for  each  of  the 
models.  The  lift  coefficient  (Cg^)  for  a  particular  angle  of 
attack  (a),  flap  setting,  and  model  was  obtained  by  adding 
the  forces  from  the  raw  front  lift  and  rear  lift  tables  and 
then  dividing  by  the  dynamic  pressure  times  the  wing  span 
times  the  chord  length  (qbc) .  Since  lift  data  for  changing 
angles  of  attack  were  obtained  for  flap  deflections  of  0“ , 

4®,  and  8®,  this  produced  three  versus  a  plots  for  each 
model.  Each  plot  intersected  the  a  axis  where  the  lift  co¬ 
efficient  was  zero,  this  determined  the  value  of  the  zero- 
lift  angle  of  attack  (ulq^  ®  specific  flap  deflection  (n) 
of  0®,  4®,  or  8®.  Three  changes  in  zero-lift  angle  of  attack 
divided  by  the  flap  deflection  (Auj^q/d)  were  then  obtained  by 
subtracting  olq  from  aLQ  then  dividing  by 

4®,  ulo  from  ulq  at  8®n  and  dividing  by  4®,  and  final¬ 

ly  OLo  at  0®n  from  olq  at  8®ri  and  dividing  by  8®.  These 
values  of  were  then  plotted  versus  the  value  of  the 

flap-chord  ratio  (E)  for  the  model  data  being  reduced  and 
against  the  theoretical  curve  from  Fig.  2.  This  gave  Fig.  11 
which  compares  experiment  to  theory  for  Aolq/h  versus  E. 

The  raw  wind-tunnel  data  was  also  reduced  to  obtain  the 
change  of  moment  coefficient  about  the  aerodynamic  center 
divided  by  flap  deflection  (AC^jj^/n)  for  each  model.  By 
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summing  moments  from  the  rear  lift  and  front  lift  forces  about 
an  assumed  aerodynamic  center  of  .245c  and  then  dividing  by 
qbc*,  the  was  obtained.  This  was  done  for  data  at  zero 

degrees  angle  of  attack  and  varying  flap  deflections  to  ob¬ 
tain  a  graph  of  Cj^ac  versus  n.  was  then  obtained  by 

simply  measuring  the  slope  of  the  C^ac  versus  n  line.  This 
value  was  then  plotted  for  the  particular  model  falp-chord 
ratio  against  the  theoretical  curve  in  Fig.  4  to  produce  Fig. 
12  which  again  compares  experiment  to  theory. 

Comparison  of  Experimental  Data  to  Theoretical  Predictions 

Figure  11  shows  the  experimental  data  plotted  with  the 
theoretical  curve  for  AoLo/n  versus  E.  Here  the  experimental 
data  compares  favorably  with  the  theoretical  curve  from  both 
a  trend  and  numerical  standpoint.  However,  for  the  small  flap 
deflections  used,  only  one  experimental  curve  was  expected  in¬ 
stead  of  three.  This  seems  to  imply  that  the  data  could  be 
biased  in  some  manner. 

Figure  12  shows  the  experimental  data  plotted  with  the 
theoretical  curve  for  versus  E.  Here,  the  experiment¬ 

al  data  indicates  the  same  general  trend  as  the  theoretical 
curve  but  the  positive  numbers  for  E  *  .6,  .7,  and  .8  are  not 
possible  unless  again  the  data  has  been  biased.  Since  it  has 
already  been  established  in  Section  II  that  the  theoretical 
curves  are,  for  all  practical  purposes,  the  same  for  cambered 
and  uncambered  flapped  airfoils,  it  is  reasonable  to  assume 
that  the  experimental  data  gathered  for  a  cambered  flapped 
airfoil  should  be  almost  the  same  as  the  uncambered  data.  The 
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Fig.  12  Experiment  versus  theory  for 

versus  E. 

mac 
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experimental  data  in  both  Figs.  11  and  12  are  overall  lower 
than  the  results  given  by  Schlichting  (Ref  4:65)  for  an  un¬ 
cambered  airfoil  with  hinged  flaps.  This  also  indicates 
that  the  data  could  be  biased  in  some  manner. 

Corrections  for  the  experimental  data  were  not  found. 
Therefore,  the  experimental  data  presented  only  serves  to 
confirm  the  general  trend  of  the  theoretical  predictions. 
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V  Conclusions  and  Recommendations 


Conclusions 

The  experimental  data  gathered  during  the  investigation 
was  biased  in  some  manner  to  produce  results  which  at  best 
can  only  be  used  to  confirm  the  trend  of  the  theoretical 
predictions  of  the  change  in  the  zero- lift  angle  of  attack 
divided  by  flap  deflection  (Aaj^Q/ri)  versus  the  flap-hinge 
location  (E)  and  the  change  in  the  moment  coefficient  about 
the  aerodynamic  center  divided  by  flap  deflection  (AC^^r/^) 
versus  E.  The  data  can  not  be  used  to  meet  the  thesis  ob¬ 
jective  of  predicting  Aolo/h  versus  E  and  AC^^r/n  versus  E 
for  a  NACA  4412  airfoil. 

While  there  is  no  direct  evidence  as  to  what  was  biasing 
the  experiment,  there  are  several  indirect  items  which  point 
to  a  major  problem.  That  major  problem  would  seem  to  be  the 
rear  lift  bar.  If  we  assvime  the  models  were  reasonable  repro¬ 
ductions  of  the  4412  airfoil,  then,  a  downward  force  was 
produced  to  counteract  the  airfoil's  lift.  The  downward  force 
had  to  be  occurring  because  from  Ref  5  the  zero- lift  angle  of 
attack  for  a  4412  airfoil  is  approximately  3.9®,  while  this 
investigation  had  a  zero- lift  angle  of  attack  of  approximately 
2®  as  demonstrated  from  the  versus  a  graphs  in  Appendix  H. 
This  is  a  large  difference  which  implies  that  something  was 
shifting  the  versus  a  graph  to  the  right /producing  a  down¬ 
ward  force.  That  much  change  could  come  from  another  lifting 
surface  which  the  rear  lift  bar  definitely  was  as  shown  by  the 
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cross  sectional  view  in  Fig.  7.  The  rear  lift  bar  was  also 
located  in  a  position  where  it  could  be  influenced  by  downwash 
changes  from  changes  in  the  amount  of  flap  deflection.  This 
could  account  for  the  large  changes  in  the  slope  of  the 
versus  a  graphs  as  the  flap  was  deflected.  Take  for  example 
the  model  with  E  =  .8,  the  slope,  in  units  1 /radian,  of  the 
versus  a  graph  goes  from  5.809  with  0“n  to  5.602  with  A®ri»  and 
to  5.204  with  8®ri.  In  other  words,  the  amount  of  flap  deflect¬ 
ion  changed  the  forces  on  the  rear  lift  bar  which  in  turn 
effected  the  model  results.  The  statements  in  this  paragraph 
imply  that  the  rear  lift  bar  was  influencing/biasing  the  exper¬ 
iment.  Unfortunately,  no  way  has  been  found  to  make  acceptable 
corrections  involving  the  rear  lift  bar  from  the  data  gathered 
during  the  investigation. 

Recommendations 

If  the  investigation  is  accomplished  again  using  the  same 
basic  airfoils  and  model,  I  would  like  to  recommend  some  changes 
which  may  prove  useful.  First,  eliminate  the  rear  lift  bar  by 
modifying  the  airfoil  to  accommodate  a  sting  for  connecting  the 
rear  lift  wire  to  the  model  or  by  replacing  the  rear  lift  bar 
with  one  of  cylindrical  shape.  Both  changes  have  inherent 
problems  but  the  rear  lift  bar  used  in  this  investigation  is 
unacceptable  and  must  be  changed  or  eliminated.  Second,  modify 
the  side  plates  to  make  them  larger  and  round.  Although  no 
evidence  exists  as  to  any  effect  the  side  plates  had  on  the  in¬ 
vestigation,  the  plates  are  aligned  Into  the  free  stream  differ¬ 
ently  when  the  angle  of  attack  Is  changing  which  could  be 


causing  different  forces  on  the  plates  themselves  and  there¬ 
fore  the  model  as  the  angle  of  attack  changes.  Third,  cut 
down  the  flap  deflection  mechanism  so  that  only  flap  deflect¬ 
ions  of  0“  to  12“  can  be  accomplished.  This  provides  more 
than  enough  data  points  to  plot  and  eliminates  wasted  effort 
since  the  tuft  studies  showed  that  little  data  is  gained 
above  12“ti  before  separation  on  the  upper  surface  occurs.  Also 
round  off  the  edges  of  the  flap  deflection  mechanism  and  cut 
down  the  size  of  the  mechanism  along  the  attachment  holes  so 
that  the  mechanism  disturbs  the  flow  as  little  as  possible. 
Fourth,  find  someway  to  seal  the  gaps  on  the  whole  airfoil  be¬ 
sides  modeling  clay.  The  clay  can  never  be  put  on  the  same  way 
twice;  tape  or  wax  might  be  two  other  possibilities.  Finally, 
the  manner  in  which  the  investigation  is  conducted  could  be 
changed  to  help  produce  successful  results.  The  investigation 
could  be  conducted  at  a  higher  dynamic  pressure  than  3  inches 
of  water  to  provide  a  Reynolds  number  of  the  same  magnitude 
that  Abbott  and  Doenhoff  had  for  their  experiments  with  the 
A412.  This  may  involve  changing  the  side  plates  to  a  different 
material  than  plexiglass  or  providing  stiffeners  for  the  plex¬ 
iglass  to  accommodate  the  increase  in  tunnel  speed.  It  would 
be  very  useful  to  find  the  drag  and  lift  of  any  side  plates 
through  wind-txjnnel  testing  so  that  they  can  be  accurately  ac- 
cotmted  for  during  the  investigation.  When  taking  data,  angle 
of  attack  changes  could  be  accomplished  in  one  degree  Increments 
and  with  flap  deflections  of  0“,  2“,  4“ ,  6“ ,  and  8“ .  This 
increase  In  data  should  result  in  more  accurate  curves  and 
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Appendix  A 


Analytical  Solution  for  ulq  and  Cmac 


Thin-airfoil  theory  gives  ulq  "  ^  (cos  6  -  1)  de  and 

^mac  “  "JyO  Hx  •  These  equations  can  be 

analytically  integrated  for  an  NACA  4412  airfoil  to  give 
numbers  for  a^o  and 

The  camber  line  of  a  4  digit  NACA  airfoil  consists  of  two 
parabolic  arcs  tangent  at  the  camber  line’s  maximum  ordinate. 
The  equations  defining  the  camber  line  are 

i=^[2pi-[ir]  <18, 

forward  of  maximum  ordinate  and 


z  _  m 
c  (l-p) ^ 


<1-  2p)+2p  I -[!]'] 


(19) 


aft  of  the  maximum  ordinate  where  m  equals  the  maximirai  ordin¬ 
ate  of  the  camber  line  expressed  as  a  fraction  of  the  chord 
and  p  equals  the  chord  wise  position  of  the  maximum  ordinate 
in  tenths  of  the  chord.  The  4  digit  NACA  airfoil  numbering 
system  gives  the  values  of  m  and  p.  The  first  digit  divided 
by  100  is  m  and  the  second  digit  divided  by  10  is  p.  (Ref  5: 


113,114) 

The  integral  equations  from  Section  I  to  be  solved  are 


and 


^mac  “  ^Jo  ^  (cos  26  -  cos  6)  d0  (13) 

The  only  other  piece  of  information  needed  to  accomplish  the 

integration  is  the  fimction  relating  x  and  6  which  is  x  * 

.5c(l  -  cos  0)  .  Thus,  finding  ^  from  equations  18  and  19 

dz 

and  substituting  for  ^  and  x  in  equation  12 

t/e  ®  <20) 

where  3  is  the  angle  corresponding  to  the  maximxjm  ordinate. 
Using  the  4412  values  for  m  and  p,  and  substituting  for  x 
again 

[.2  -  .25(1  -  cos  0)]  (cos  0  -  1)  d0 

+  /e  ](cos  6  -1)  doj  (21) 

The  angle  B  corresponding  to  m  and  p  can  be  found  easily  by 
substituting  in  the  value  for  x  at  the  maximum  ordinate  into 
the  equation  x  **  .5c(l-cos0)  .  Thus,  for  the  4412,  .4c  = 

.5c(l-cosB)  which  implies  that  cos  6  =  .2  and  therefore 
B  *  1.3694384  radians.  By  expanding  the  integral  and  putting 
the  numerical  value  of  B  into  equation  21 


Lj 
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«• 


*L0 


r  ^  A .  oo 

4/0  [• 


1.369 

2(cos  0  -  1)  -  .25<2cos  0  -  1  -  cos^  0)]  d0 


II 

[-TT  (cos  0  -  1)  -  4  (2  COS  0  -  1  -  COS^  0) 


1.369 


After  integrating  equation  22 


)]d6] 


(22) 


«L0 


■-[[ 


2  (sin  0  -  0)  -  .5  sin  0  +  .250  +  + 


+ 1  ( s  in  0  -  0)  -  •^sin0+^  +  '^ 


[|.sin^]]|’ aej  (23) 


By  putting  in  the  limits  for  0 ,  we  arrive  at 


“LO 


-1 


(.227795)  =  -.072509A  radians  =  -4.1544826  degrees 


Following  the  same  steps  as  for  ulq 
1.369 


if/* 

g  =  yI /q  [.  2(cos  20  -  cos  0)  -  .25(1  -  cos  0)  (cos  20  -  cos  0)]  d0 

(cos  20  -  cos  0)  - -5  (1  -  cos  0)  (cos  20  -  cos  0)j  (24) 

*'1.369  J 

c  ■  -  -zsf 


0  ,  sin  20 


5 - sin0+ 


sin  20 


.125[! 


sin  0  + 


sin  30 


-[ifE 


-  sin  0] 


sin  20 


-it 


sin^  _sine+  1+  ]  +.0555556  [sin  0 


Cmac  -  ^(-.2124781)  -  -.1062391 


(25) 
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This  appendix  contains  the  listing  of  the  computer  program 
used  to  find  the  theoretical  predictions  for  zero- lift  angle 
of  attack  (olq^  »  moment  coefficient  about  the  aerodynamic  cen¬ 
ter  (Cjjjac)  »  change  in  aLO  divided  by  flap  deflection  (AoLQ/n) 
versus  the  flap-chord  ratio  (E) ,  and  versus  E  as  given 

in  Section  II.  Comment  cards  in  the  program  define  the  vari¬ 
ables,  parameters,  and  limitations.  The  program  used  n\americal 
integration  to  calculate  the  predictions. 

The  program  did  not  attempt  to  keep  computer  time  to  a 
minimum  or  provide  a  program  which  needed  a  minimal  number  of 
segments  to  obtain  accurate  results.  The  objective  was  simply 
to  obtain  accurate  results  as  demonstrated  by  the  comparison 
between  the  analytic  solution  and  the  program  solution  for  aLO 
and  at  the  top  of  page  6.  It  should  be  pointed  out  that 

the  method  used  in  this  program,  the  trapezoidal  rule  of  inte¬ 
gration,  does  have  inherent  errors  and  requires  a  relatively 
large  number  of  segments  to  achieve  accurate  results.  The 
accuracy  could  be  Improved  and  the  number  of  segments  needed 
to  compute  an  accurate  solution  could  be  reduced  by  using  a 
different  ntunerlcal  integration  method.  Since  the  mean  camber 
line  of  the  NACA  A412  airfoil  is  a  parabolic  arc,  Simpson’s 
Rule  which  uses  parabolic  arcs  could  be  used  to  produce  better 
accuracy  with  fewer  segments  when  compared  to  the  results  of 
this  program. 


3  NOV  e?  POli^tOO  HirfPlS  FORTRAN  77  SA'J  OoTT^TZING  cc^^pilfr 
MODULE  NAME:  *“AIN* 


01 ai tt-on 


1 

2 

3 

a 

5 

6 
7 
6 
9 

10 

11 

12 

13 

la 

15 

16 

17 

18 
IR 
20 
21 
22 
23 
2a 

25 

26 

27 

28 
20 

30 

31 

32 

33 
3a 

35 

36 

37 

38 

39 

ao 

01 

02 

03 

00 

05 

06 

07 

08 

09 

50 


PROGRAM  FLOP 

DESIGNED  BT  HOWARD  P«»1CE»  »UG  02,  THFSTS  ^PPK  0'  VORVi-c  hTwGEP 
FLAP. 

Integer  i, j,K,L»LE.LfTA,N,NA,NAc»oo,NE,NFs,NP,NPTs.LL 

real  AF,  AF\,  ALO,  AUQ2,  A3.Aa,CMAC,r»-Ar2,neLl  ,DFL2»0ELlT,t'FL2T,»'S,Dy, 
$OXA’,OXT,n2,D70X,0ZN,OZT,E,FT,FTA,FC,FCl  ,FC?,>'/P/PX,P7,SA,SA.-i,SF, 
jSL»STH,S2TH,THETA,X,yc,  XCn,xf,xn,X1 ,y2,XX,7,7N,Z? 
dimensions  FOR  PCINT/SFC‘'£*'T  PROPERTIES  OF  ATRFOTL  CA'PEPlInF  ANO 
0  OF  FLAP  CHORD  RATIOS.  CHAt.'GE  OImE'‘  SlPf'S  FCP  MORE  I  '*-*'  «00 
SEGMENTS  OR  100  FLAP  CHORD  HATIOS. 

DIMENSION  AF(80n,AFN(A01),tFL1(lOl.lO),r>EL2(10l,l'').OS(»iO),D7(«0 
Sl).07PX(R00),E(10n,LETA(ie),P*(«01  ),PZf  son.  theta  (PA  11,  SA  {8011,50 
sN(8oi)  tSF(8oi)(SL(8oii«y(8on,y'«(8on.z(noi )  ,ZN(8on 
INPUT  «(NA)  OF  AIRFOIL  C A.^PEPLInES  TP  PE  AMAI  TZ£n 
REA0(7»a)NA 

MAIN  control  loop 

DO  16  Kst,NA 

INPUT  0  OF  SEGMENTS(N),  AIRFOIL  *(hACA001,  ChOpO  UFjCTM{KC) 
REAO(7,*)H,NACAaO«XC 

*  OF  POINTS(NPTS),CAMPER(M),  MAX  C*mREP  L OC ATIp; (P ) 

NPISsN+l 

GET  camber  from  nACAOO  • 

MsNACAOO/IOOO. 

MsAINT(«) 

Ma.OlPM 

GET  max  camber  LOCATlCN  FRO“  NACAUt)  « 

PsNACAOO/lOO. 

PsAINT(P) 

Ps,|*(P-(M«lOOO.)) 

POTNT(NP)  OF  max  camper 
NPS(P*N)>1 

X  COMPONENTIOXI  of  every  SFGMfHT 
OXsxC/N 

NACA  0  DIGIT  AIRFOILS  CONSIST  OF  2  PARABOLIC  ARCS  TANGF».T  AT  ^AX 
CA“8ER  point.  00  LOOP  1  CQmpuTFS  »  p  TC  ^P  ':51.;G  1  »'«C  •■^rlLF  pO 
LOOP  2  COMPUTES  FRO*  nP  TO  TPAILT-G  EDGE  FpI»'T('P(S)  URI*G  2mo 
ARC.  COORDINATE  SYSTEM  MAS  LEAPING  FOSE  “OlftT  aT  ORIGIN  ANO  MPTS 
ON  X  AXIS  AT  XC  FRO"  CRIGIM. 

PO  1  J=1»NP 

COOPOINATF(X(J1)  or  POINT  ON  x  AXIS.  XfPXlJl)  X(.!)  T5  OF  XC. 

x(j)s(j-t)«nx 

pxf J)=x<jj/xr 

IFtP.GT.O.lTHEM 

FOR  cambered  AIRFOTL,  ARC  FCUATIO'  FT?  PS  KFiGr  T(7(  J.i)  op  ' 
POINT  On  CAnBERLInE  I?  X(P7(J))  OF  XC. 

P2(JlsM*(  (2. *P*PX(jn-(PXl. ’)*•?) 

ELSE 

symmetrical  AIRFOIL. HEIGHT  OF  rA»/JjFPLIN1.  IS  o. 

PZ(J1»0. 

FNO  TF 
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3  N'OV  b?  ?0;19:00  Ht«RlS  FOPTPaW  77  SAU  nPTlHIZThr,  CC“PTLFR 
ROOULE  name;  FUaP 


Ota!<i-On 


51: 
52: 
53: 
5p: 
55: 
56; 
57; 
58; 
5R; 
60: 
61: 
62; 
63: 
66: 
65; 
66: 
67: 
68: 
69: 
70: 
71: 
72: 
73: 
7a: 
75: 
76: 
77: 
78: 
79: 
80: 
8l: 
82: 
83: 
84: 
85: 
86: 
87: 
88: 
89: 
90: 
91 : 
92: 
93: 
94: 
95: 
96: 
97: 
98: 
99: 
100: 


C 

C 

C 

c 

c 


c 


c 


c 

c 

c 

c 

c 

c 


2 


Z(J)=XC«P7{J) 

PlSTANCE(SL(JnFRf)M  POINT  t  TO  POINT  J, 
SL(J)sSQRT((X(J)**2)*(Z(J)**p)  ) 

4NGLF(SA(J))  .vITH  X  AXTS  OF  IINE  FRO«  POir-T  1  TO  POJVT  J 
SA  IN  RADIANS 

CHANGE  IN  2(07(J1)  FROM  POTNT  TO  POINT 
IFfJ.EO.DTHEN 
FRONT  CONDITIONS 
SAf J)sO. 

OZIJlsO. 

ELSE 

SA(J)=ATA:4(Z(JI/Xf  J)1 
OZ(J)aZ(J)-Z{J-l) 

END  IF 

X(J)  TRANSFORMATION  ANGLEITHETA ( J) )  IN  PlCIANS 
THETA  (J)=AC0Sn,-(2.*XfJ)/XC>) 

CONTINUE 

00  2  JsNPfltN 

X(J)exU-l)t0X 

PX(J)=X(J)/XC 

PZ(J)=P»(l,-(2.*P)+I2.*P*PX(J))-(PX(J)*A2n/f (« .-P)»»2) 
2(J)»XC»PZ<J) 

0Z(J)s2(J)-Z(J-n 

THETA(J)aAC0S(1.-(2.*X(J)/XC)) 

SLU)aSORT((X(J)**2>  +  (Z(J)**37> 

SA(J)=ATAN(Z(J)/X(J)) 

continue 
rear  conditions 

X(NPT3)sXC 

THETA(NPTS)=3.14159Z65 

SU(M>TS)=XC 

SA(NPTS)=0. 

ZINpTSJsO. 

OZ{NPTS)sZ{NPTS)-Z(N} 

INITIALIZE  ANGLE  OF  7ER0  LTFT(ALn)  AnP  C^EFFICIEM  OF  -PHEMT 
ABOUT  AEROOTNAPIC  CENTFRICmaC) 

ALOsO. 

CNAC=0. 

NUMERICAL  SUMMATION  PNOCESS/IOOP  FpR  ALO  Ar'D  C^^AC 
00  3  I31»N 

SLOPE(CZCX(I))  OF  EACH  SECMEMT 

Dzox(i)=oz(i+i)/nx 

LENGTHCOSdIlOF  EaCN  SFGmenT 
0S(I)=S;RT(r02(I+l)**2TtICX**2l) 

intermediates  in  ALC  AFP  CMAC  FfiiATlCNS 

3TH«SIN(THETA(T4m«aiMThrTA(T)  ) 

S2ThsSIn(2.*TH«:TA(U1  ))-SIV{?.»TmETa(T>) 
ALOs(OZOxn)*(STHfTHFTA(n-TMET:  CI  +  nil  +  ALO 
cn4Cs(0zox(I)«(-sth+ts?tm/?.)) i+rwAt 
CONTINUE 


01  4l<i-0n 


t. 


3  NOV  82  20:10:00  HAfiOjS  FOOIPan  77  SAu  OPTt»'tZlN»;  CO^PILFS 
MODULE  name:  flap 


101:  C 

102: 

103:  C 

104: 

105:  C 

10b: 

107: 

108: 

109: 

no: 
111:  C 
112:  C 
113: 
114: 
115: 
116: 
117:  C 

lie: 

119: 
120: 
I2i:  C 
122:  C 
123: 
124: 
125: 
12b: 
127: 
128:  c 

129:  C 

130: 

131: 

132: 

133:  C 

134i 

135: 

136:  C 

137:  C 

138: 

139: 

140: 

141: 

142: 

143: 

144: 

145: 

146:  C 

147: 

148: 

149: 

150: 


final  ALO  in  cegpees 

ALO=(-ALO/3.141592651»(ieO,/3. 14159265) 

FINAL  CMA.C 
CmACsCMaC/2. 

OUTPUT  KITHOUT  FLAP  DEFLECTION 
PRINT* ("I")  ' 

PRINT* CO", 33X,"ALPHA(L0)  FOUALS  ",F10.7,"  .  CMAC  EQUALS  ", 

$  FIO. 7, ".’)*, ALO, CMAC 

PRInT*(38X."  fixed  FLAP  CHORD  RATIOS-)' 

PRINT*(38X,"  varying  flap  DEFLECT  IONS' ) * 

L  =  dimension  variable,  LE  =  LOOP  COUNTER  TO  ChA*'GF  ^LAP  CHORO 
RATIO  E, 

L=0 

DO  13  LEs0,100 
LsL+1 

E(L)s.01"LE 

NE  =  HINGE  POINT,  NFS  =  •  OF  FLAP  SEGPFNTS. 
NEsNPTS-((LE*N)/100) 

NF8=NPTS-NE 
DO  6  1=1, NFS 

OXT  s  X  PROJECTION  OF  1  FL*P  SEGMENTS, 

OZT  •  Z  PROJECTION  OF  I  FLAP  SEGMENTS, 

OXTsIaOX 

OZTsO, 

00  5  J81,I 

OZTaDZ(J»NE)+OZT 

5  CONTINUE 

SF  s  OIST  FROM  HINGE  POINT  TO  A  FLAP  POINT,  A  FLAP  LINE, 

AF  s  angle  NTTH  Z  axis  SF  “AXES, 

SF(I)=SORTf {DXT«*2)t(07T**2)) 
AF(I)sI,5707963*ASIN(OZT/SF(I)) 

6  CONTINUE 

XE  =  FLAP  CHORD  LENGTH  X  PROJECTION,  FC  =  FLAP  CMQRD  LENGTH, 
XE3E(L)aXC 

FC=SQRT((XE«*2)+{2(NE)**2)) 

NFS  s  END  POINT  CONDITIONS,  ^E  =  1  means  ENTIRE  AIRFOIL  IS 
A  FLAP,  ME  *  NPTS  means  the  AIRFOIL  HAS -NO  FLAP, 

SF(NFS)aFC 
IF(NE,Ea,I)THEN 
AF(NFS)sl, 5707963 
ELSE  IF(NE,EQ,NPTS)THEN 
AF(NFS)st,57079i>3 
ELSE 

AF (NFS) a  1, 570  7963- AT AN(Z(NE)/XF) 

END  IF 

SETUP  OUTPUT  Table, 

PRINT*  ("0",20X, "FOR  A  FLAP  C“i>RC  RATIO  OF  -.Fb.S, 

S  ",  THE  FLAP  REFLECTS  ARO'iT  R?iI'’T-,I0, 

S  ■  AND  GIVES  the  follCh^tno  Data : • )  *  ,e  (L ) , nE 

PRINT* (22X, "ETA", 2X,-ALPHA(L0)", MX, "rPAC'.SX, 
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3  NOV  82  20:i‘>:00  MABHIS  F0»T»AN  77  SaiJ  ObJ  I«T Z T«r.  (“OWPILES  Q\U\H-0n 

MODULE  na«E:  flap 


JSIS  s 

152: 

153: 

iSa: 

15S:  C 
l5fc:  C 
157: 

158: 

IS*!:  C 
160:  C 
161:  C 
162:  C 
1631  C 
1641 
165: 

166: 

167: 

168: 

16‘?:  C 
170: 

171:  c 

172: 

175:  C 
17<»:  C 
175:  C 
176: 

177: 

178: 

170:  c 

180:  c 

181:  C 
182: 

183: 

tea: 

185: 

18(,: 

187: 

188: 

189; 

190: 

191 : 

192: 

193: 

190: 

195:  7 

196: 

197: 

198: 

199: 

200:  8 


•DEL(ALPHA(L0) VeTA«,2»/«0ELfC“Ar)/ETA') • 

LL-0 

DO  12  JJ=1.10 
LL=LL4-1 

LL  =  OIMENSION  VARIA9LF,  JJ  sLPOP  CONTROL  FOB  iy.nu^'T  OF  FXAP 
deflection.  LETa  =  DEFLECTION  TN  CEGBEES,  ETA  =  BAOIAfjS. 
LETA(LL)=J.I 

ETA=LETA(LI.)*3.10159?6';/18!'. 

GEOMETBY  FINOS  NE«  COOBOfiATFS  FOB  FLAPPFD  AIRFOIL. 

8£T«eEM  flap  CHOBO  AMO  X  *xIS»  Xt  s  o^OJEC TI 0»  ONTO  X  AXIS 
AFN  =  NEV<  ANOLF.  a  flap  line  nAKES  ^ITH  Z  AXIS,  A3  s  angle 
OF  FLAP  CHORD  ABOVE  X  AXIS.  FCl  =  LEnGTH  CF  e’LAP  CHOBO 
ABOVE  X  AXIS.  FC2  *  LENGTH  OF  FLAP  CHORD  0ELO.1  X  AXIS, 
AFN(NFS)sAF(NF8)-ETA 
A3s1.5707963-AFF!(NFS) 

X1sZ(NE)/TAN(A3) 

FC1=S0RT( (X1**2)+(Z{nE)**2I) 

FC2=FC-FCl 

22  a  PROJECTION  ONTO  Z  AXIS  CF  FLAo  CHORD  BELOH  <  AXIS. 
72aFC2*3IN(A3) 

X2  s  PROJECTION  ONTO  X  AXIS  OF  FLAP  CHORD  RELOY  X  AXIS 
X2aZ2/TAN(A3) 

X3  a  PROJECTION  OF  NEW  CHORD  LINE  ONTO  X  AXIS.  XCN  a  LENGTH 
OF  NEW  CHORD  LIME  AFTER  DEFLECTION.  AO  a  ANGLE  BEIWEFn  NE.. 
CHORD  LIME  ANO  X  AXIS. 

X3aXltX2'»X(NE) 

XCNaS3RT((X3**2)*{Z2*«2)) 

Aa3ATAN(Z2/X3) 

COMPUTE  NEW  COOROINATES  AND  angles  with  TFF  NE'  chord  line 
being  the  new  X  AXIS,  LOOP  7. TAKES  CARF  OF  POINTS  I'p  TO 
hinge  point, B  and  9  DO  POINTS  FROM  HImGE  PT  TO  E' D  POINT. 

DO  7  Jrl.wE 

IFIJ.EO.llTHFf!. 

SAN( J)aO. 

XNdlaO. 

ZNClIaO. 

ELSE 

SAMf.t)r9A(  J)»A4 

X-I<  J)rSL(J)*COS(SAN(J)> 

Z;.(J13SL(J1*STM(.SAN(J)) 

OXMsXNI  J)-XM(  j-n 
OZNaZNI  JI-7MfJ-n 
OZOXlJ-naOZN/OXN 
END  IF 
CONTINUE 
00  B  J81,NFS 

AFN(  J)8AF{  J)4A«-F.7A 

XN(JfNE)aXN(NEUfSF(J)*STN(»FN(J))) 
rN(J*NE)aZN{Nf )-ICnS(AFN(J))*sr(J)) 

CONTINUE 
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XNCNPTS)=J(CN 
2,V(NPTS)=0. 
ro  0  J=MFtl /UPtS 
nx6sXN(  J)-XM(  J>n 
0ZNs7\'(  Jl-2w(j-.n 
rzox(j-nsn2^'/oxf! 

0  COf^TlN'JE 

N£/r  Tf^A^(SF0R^•ATTf1.M  AN6I.F?>. 

CO  10  J=1,F' 

ThETA(  JlsAfOSCI  .-(2.*XA>(J)/XCN>) 

10  CONTIW'JE 
AU02=0. 

ChAC2=0. 

ThETA(NPTS)=3. 10150265 
NEW  ALO  ANO  CKAC. 

00  11  1=1,6 

STHsSlN{THrTA(I+1))-SlW(ThFTA(I)) 

S2THsSIN(2.*THETA(I  +  m-SlM?.*THETAf  ID 

AL02a(0Z0X7I)*(STH+TKETACl)-THETA(ltl)>)»ALO2 

CK4C2=<0Z0X(n<*(-STH»CS2TH/2.)))»CPAC2 

11  continue 

ACCOUNT  FOB  NEW  ALO  BEING  fP“PMTFO  IN  OIFFEMEM  COORDINATE 
SYSTEM  AND  CPNVfOT  TP  0E6BEES. 
AL02a((-AL02/3.1fll50265)-A«i)«(iBO./3.14150265) 

CMAC2*rMAC2/2. 

FIND  CHANGE  IN  ALO/FLAP  DEFIFCTIDK  FOR  E. 
DEl.l(L,LL)*(AL02-AL0)/lETAaL) 

FIND  CHANGE  IN  CMACZFLAP  DEFLECTION  FOR  E. 

0EL2(L,Ll  )»(r.MAC2-C«AC)/ETA 
OUTPUT 

PRINT* (22X, I2,3X,F9.4,3X,F7.4,4X,F9.4,  U>X,F9,4i ' ,LFT*(LL), 

S  ALn2,CMAC2,0FLl (L,LL),CEL2fL,LL) 

12  CONTINUE 

13  CONTINUE 

OUTPUT  IN  A  OIFFERFNT  FORM, 

PRINT* ("1") • 

PRINT*  (SOX, **  FIXED  FLAP  PCFLFCTinN")  * 

PRINT* Oax,"  VARYING  FLAP-CHORD  RATIO") • 

DO  15  1=1, LL 
PRINT* ("0")  • 

PRINT* (22X, "FOR  A  FLAP  DEFLECTION  OF  ", lo, "DEGREES i")' ,LETA(  I) 
PRINT* (22X"E(FLAP-CHn»ft  RATIO)", ?X, "DEL  ALPHA(LC)/FI A",3X, 

%  -OF.  CNAC/ETA")* 

DO  14  Js),L 

PRINT*  (30X,F5.3,10X,F6.3,12X,F5.3)  *,F(.l),l5rLl  (.i,T),rEL2(J,I) 

14  CONTINUE 

15  CONTINUE 

16  CONTINUE 
ENO 
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Appendix  C 

Model  Geometry 

The  side  plates  were  two  and  one-half  feet  square  with 
each  corner  rounded  to  an  arc  with  a  four  inch  radius.  The 
rear  lift  bar  dimensions  were  1.15  inches  in  chord  length 
and  .5  inches  in  thickness,  see  Fig.  7.  See  Abbott  and  Von 
Doenhoff  (Ref  5:113,11A,  and  41A)  for  cross  section  dimen¬ 
sions  of  the  NACA  4412  airfoil. 

Table  II  gives  the  airfoil  measurements  needed  to 
compute  Cj^s  and  for  each  model. 


TABLE  II 

Airfoil  measurements 


Model 

E 

Chord  length 
(inches) 

Wing  span 
(feet) 

Distance  from  trunnion 
screw  to  leading  edge 
of  airfoil  (inches) 

.2 

8.9 

1.992 

.3 

.3 

8.9 

1.992 

.31 

.4 

8.88 

1.992 

.31 

.5 

8.88 

1.987 

.33 

.6 

8.89 

1.992 

.3 

.7 

8.88 

1.99 

.3 

.8 

8.9 

1.99 

.3 

The  model's  angle  of  attack  was  set  by  moving  the  rear 
lift  wire  forward  or  backward  and  Increasing  or  decreasing 
the  length  of  the  rear  wire  from  the  scale  to  the  rear  lift 
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bar  attachment  fitting.  The  angle  of  balance  in  Fig.  13  de¬ 
termined  the  distance  forward  or  backward  the  rear  lift  wire 
had  to  move  for  a  specific  angle  of  attack  while  the  shorten¬ 
ing  or  lengthening  of  the  wire  was  determined  during  calibra¬ 
tion.  Table  III  gives  the  balance  angles  and  model  measure¬ 
ments,  shown  in  Fig.  13,  that  determined  them.  Table  III 
shows  that  the  balance  angles  and  model  measurements  varied 
barely  at  all.  For  that  reason,  the  first  model  tested  with 
E  =  .4  determined  the  angle  of  balance,  and  therefore  the 
distances  the  rear  lift  wire  moved,  for  the  entire  investiga¬ 
tion. 

TABLE  III 

Balance  angles  and  associated  measurements 


Model 

E 

Distances  in  Fig.  13 
(inches) 

Balance  angle 
(degrees) 

e 

f 

g=Ve*+f2 

.2 

15.737 

9.545 

18.405 

31.238 

.3 

15.776 

9.637 

18.487 

31.419 

.4 

15,71 

9.516 

18.367 

31.204 

.5 

15.726 

9.577 

18.413 

31.341 

.6 

15.684 

9,599 

18.388 

31.468 

.7 

15.737 

9.572 

18.419 

31.31 

.8 

15.741 

9.694 

18.487 

31.627 

The  distances  the  rear  wire  moved  are  given  in  Table  IV. 
The  distance  and  direction  of  movement  for  a  specific  angle 
of  attack  in  Table  IV  was  determined  by  calculating  a  new 
balance  angle  for  each  angle  of  attack.  The  new  balance 
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Fig.  13  Balance  angle  and  associated  measurements. 

e  was  the  horizontal  distance  and  f  was  the 
vertical  distance  from  the  front  trunnion 
to  rear  lift  wire  attachment  fitting  with 

g  =  \/e*+f  *  . 

angle  gave  a  new  distance,  e  in  Fig.  13,  to  the  rear  fitting 
from  the  front  trunnion.  The  difference  between  the  distance 
e  at  O'*  angle  of  attack  and  e  at  a  specific  angle  of  attack 
gave  the  distance  and  direction  the  rear  lift  wire  was  moved. 
These  distances  were  put  on  a  bar  graph  and  used  to  set  the 
rear  lift  wire  from  the  wind-tunnel  operator's  station  during 
testing  that  required  changing  the  angle  of  attack.  As  stated 
previously  the  distances  of  movement  in  Table  IV  were  used 
for  the  entire  investigation  and  came  from  the  model  with  E  ■ 
.4.  The  same  distances  were  calculated  for  all  the  models 
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and  for  all  practical  purposes  do  not  vary.  For  example,  the 
model  with  E  -  .B  gave  a  maximum  difference  from  the  distances 
in  Table  IV  at  the  angle  of  attack  of  10“  and  amounted  to  only 
-.031  inches. 


TABLE  IV 

Rear  lift  wire  movement 


Angle 

of 

attack 

(degrees) 

Balance  angle 
(in  degrees) 

Horizontal 
distance  e 
(inches) 

Rear  lift  wire  movement 

+  towards  airfoil 
away  from  airfoil 
(inches) 

0 

31.204 

15.71 

0 

-6 

37.204 

14.629 

1.081 

-4 

35.204 

15.008 

.703 

-2 

33.204 

15.368 

.342 

2 

29.204 

■39 

-.323 

4 

27.204 

16.336 

-.626 

6 

25.204 

16.619 

-.909 

8 

23.204 

16.882 

-1.172 

10 

21.204 

17.124 

-1.414 
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Appendix  D 

Flap  Deflection  Settings 

Flap  deflections  were  accomplished  by  aligning  one  of  the 
holes  (1-7)  of  the  brass  plate/flap  deflection  mechanism  in 
Fig.  6  with  one  of  the  holes  in  the  side  plate  in  Fig.  lA. 

Each  set  of  vertical  holes  were  used  with  a  particular  model 
as  shown  in  Fig.  8.  Additionally,  each  set  was  drilled  in  a 
circular  arc  about  the  flap-hinge  point  of  the  particular 
model  the  set  was  used  with.  This  allowed  each  hole  of  a  set 
to  be  spaced  at  6“  intervals.  The  holes  (1-7)  of  the  brass 
plate  were  drilled  in  a  circular  arc  about  the  hinge  point  A, 
see  Fig.  6,  of  the  brass  plate  but  were  separated  by  5®  in¬ 
stead  of  6® .  The  differences  in  the  increment  between  holes 
in  the  brass  plate,  5®,  and  a  set  of  vertical  holes  in  the 
side  plate,  6®,  allowed  flap  deflections  of  0®  to  20®  either 
up  or  down.  By  numbering  each  set  of  vertical  holes  as  shown 
in  Fig.  lA  and  using  the  numbering  system  in  Fig.  6  for  the 
holes  in  the  brass  plate,  Table  V  was  developed  to  give  the 
flap  deflection  for  particular  combinations /alignments  of 
holes  In  the  brass  and  side  plates.  Take  for  example  the  flap 
deflection  shown  in  Fig.  lA.  The  set  of  vertical  holes  being 
used  tells  you  that  the  particular  model  had  a  flap-hinge 
location  (E)  of  .A.  The  side  plate  hole  number  used  was  6  and 
the  hole  number  for  the  brass  plate  was  5.  For  no  deflection 
at  all,  both  ntsnbers  must  be  A.  So,  if  the  brass  plate  hole 
number  4  was  aligned  and  set  with  hole  number  6  of  the  side 
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plate,  the  deflection  would  be  12*  downward  since  the  side 
plate  holes  were  spaced  at  6®.  But,  since  the  brass  plate 
hole  was  5,  that  had  the  effect  of  bringing  the  flap  back 
upward  by  5®  to  give  a  flap  deflection  of  7®  downward.  The 
thinking  process  in  the  above  example  can  be  accomplished  to 
arrive  at  each  of  the  flap  deflections  in  Table  V. 


1 

2 

3 

4 

5 

6 
7 


Fig.  14  Side  view  of  flap  deflection  mechanism  in 
operation.  The  number  above  each  set  of 
7  vertical  holes  designates  the  model's 
flap-hinge  location  used  with  that  set  of 
holes.  The  vertical  nianbering  system 
designates  the  side  plate  hole  numbers. 
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Appendix  E 

Tuft  Studies 

Tuft  studies  were  accomplished  by  placing  nylon  strings 
on  the  upper  surface  of  the  airfoil,  see  Fig.  15.  Table  VI 
gives  the  results  of  the  studies  by  indicating  the  combina¬ 
tion  of  flap  setting  and  angle  of  attack  at  which  separation 
was  first  noticed.  Figures  16  through  21  determined  the  data 
for  model  with  E  =  .8  in  Table  VI.  The  figures  were  obtain¬ 
ed  by  taking  pictures  through  a  window  on  top  of  the  wind 
tunnel. 


TABLE  VI 

Angles  of  separation  from  tuft  studies 


Model 
flap -chord 
ratio 
<E) 

Angle  of  attack 
fixed  at  0* 

Separation 
angle  of  attack 
for  flap  fixed  at 

Separation 
flap  angle 

rai 

mm 

8" 

.2 

13" 

10" 

8" 

8" 

.4 

11" 

- 

10" 

6" 

.5 

12" 

10" 

8" 

6" 

.6 

O 

o 

10" 

8" 

6" 

.8 

15" 

10" 

8" 

6" 

5 


Fig.  16  8*  angle  of  attack  and  0*  flap  deflection. 

Use  this  figure  for  comparison  purposes 
since  there  was  no  separation. 


52 


Fig.  19  8®  angle  of  attack  and  A®  flap  deflection. 
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Raw  Tunnel  Data 


The  following  tables  contain  the  data  taken  directly  from 
the  wind-tunnel  tapes  for  all  seven  models.  The  tables  for  a 
particular  model  are  indexed  into  three  parts  with  part  a 
containing  the  front  lift  data,  part  b  the  rear  lift  data, 
and  part  c  the  drag  data.  The  drag  tables  contain  the  tem¬ 
peratures  in  degrees  Fahrenheit  and  the  pressures  in  inches 
of  mercury  that  occurred  during  the  testing  of  each  model. 
Angles  of  attack  (a)  and  flap  deflections  (n)  are  in  degrees 
while  the  forces  and  wind-on/static  readings  are  in  pounds. 
The  force  colximn  in  each  table  was  computed  by  subtracting 
the  wind-on  average  from  the  static  average.  The  forces 
were  used  as  explained  in  Appendix  G  to  calculate  the  exper¬ 
imental  data  points. 


Rear  Lift  Data  for  Model  with  E 


*»g.  „  a  Static  Yind-on 

Press.  Eefore-Statlc  Readings  '.Yind-on  Readings  After-Static  Readings  Average  Average 


TABLE  Xlllb 

Raw  Rear  Lift  Data  for  Model  with  E 


O  iH 
M  <0 
M  "O 
H-l  O 

X  s 

►J  o 
CQ 
< 

Eh  CO 
-*-> 
<0 
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Appendix  G 


Data  Reduction  Calculations 

This  appendix  explains  in  detail  and  by  example  how  the 
experimental  data  in  Appendix  F  was  reduced  to  obtain  the 
experimental  points  for  the  plots  of  experiment  versus  theory 
for  the  change  in  zero- lift  angle  of  attack  divided  by  flap 
deflection  (AaLo/il)  versus  flap-chord  ratio  (E)  and  the  change 
in  moment  coefficient  about  the  aerodynamic  center  divided  by 
flap  deflection  ( AC^^^/n)  versus  E,  Figs.  22  thru  28  in 
Appendix  H.  The  model  used  as  the  example  is  E  =  .7.  The  re¬ 
duction  process  for  the  AUj^q/ti  versus  E  experimental  points 
will  be  described  first. 

The  total  lift  generated  by  the  model  at  a  particular 
angle  of  attack  and  flap  setting  was  obtained  by  adding  the 
lift  force  from  the  raw  front  lift  data  to  the  lift  force  from 
the  raw  rear  lift  data.  The  lift  was  then  divided  by  the 
dynamic  pressure  times  the  wing  span  times  the  chord  length 
(qbc)  to  produce  the  lift  coefficient,  Cj,.  The  d3mamic  pres¬ 
sure  was  converted  from  3  inches  of  water  to  lbs /ft*  by 
finding  the  weight  of  one  cubic  foot  of  water  at  the  testing 
temp,  dividing  by  12  to  give  the  weight  of  one  inch  of  water 
on  one  square  foot,  and  then  multiplying  by  3  to  give  the 
pressure  of  3  inches  of  water  in  Ibs/ft*.  At  this  point, 
should  have  been  corrected  for  any  errors  induced  by  the  wind 
tunnel.  However,  the  corrections  from  Pope  (Ref  6:305-318) 


for  the  two  dimensional  wind  tunnel  were  so  small  that  they 
were  considered  negligible  for  these  models.  was  then 
plotted  versus  its  angle  of  attack  for  a  particular  flap  de¬ 
flection.  This  process  was  accomplished  for  the  three  flap 
settings  of  0®,  4®,  and  8®  at  angles  of  attack  -6®,  -4®,  -2®, 
0®,  2®,  4®,  6®,  8®,  and  10®.  This  gave  three  plots  for 
versus  a  for  a  model  with  a  particular  E.  Each  plot  had  an 
angle  of  zero- lift  and  by  subtracting  them  appropriately  the 
changes  in  the  angle  of  zero-lift  (Aolq)  were  found.  Divid¬ 
ing  the  Aolo®  appropriate  change  in  flap  deflection 

then  gave  three  Ao^^q/ti  experimental  points  for  a  model  with 
a  particular  E.  Table  XIV  gives  the  numerical  values  for  L, 
qbc,  and  for  the  model  with  E  =  .7.  Figures  22a  through 
28b  give  the  values  for  olq  seven  models. 

ACniac/’l  versus  E  experimental  points  were  obtained  in  a 
manner  similar  to  Aulq/d  versus  E  points.  The  calculations 
about  to  be  presented  were  initially  done  with  an  assumed 
aerodynamic  center  of  .25c.  This  resulted  in  a  positive  re¬ 
sult  for  ACjjjgg/n  for  E  =  .7  and  was  the  first  indication  that 
the  data  was  biased  in  some  manner.  Calculations  for  the 
aerodynamic  center  from  the  data  at  0®  flap  deflection  gave 
an  location  of  approximately  4  inches  which  was  not  poss¬ 
ible  for  the  9  inch  4412  airfoil.  In  fact,  for  the  low 
Reynolds  number  of  the  airfoil  in  this  experiment,  it  was  more 
likely  that  the  aerodynamic  center  was  in  front  of  the  quarter 
chord.  Since  no  reliable  data  for  the  Reynolds  number  of 
approx  5x10’  could  be  found,  an  assumed  location  of  .245c  was 
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TABLE  XIV 

Calculations  Tor  Model  with  E  =  ,7 
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t. 


made  from  trends  of  the  aerodynamic  center  location  in  Abbott 
and  Doenhoff  (Ref  5:A89).  The  calculations  were  then  repeated 
with  .245c.  The  forces  for  the  moment  were  the  rear  lift  and 
front  lift  with  the  front  lift  producing  a  positive  moment  and 
the  rear  lift  producing  a  negative  moment  when  the  rear  lift 
was  positive.  The  moment  arm  for  the  front  lift  was  found  by 
taking  .245c  of  the  chord  length  plus  the  distance  from  the 
trunnion  screw  to  the  model,  2.476  inches  for  E  =  .7.  The  rear 
moment  arm  was  found  by  subtracting  the  front  lift  moment  arm 
from  the  horizontal  distance  e  as  given  in  Table  III,  thus  the 
rear  lift  moment  arm  for  the  E  =  .7  model  was  13.261  inches. 
After  finding  the  moment  for  a  particular  flap  setting  at  0“ 
angle  of  attack,  it  was  then  divided  by  qbc*  to  give  C^ac* 

Table  XIV  gives  the  calculations  for  with  E  =  .7.  This 

value  of  C^jjac  then  plotted  against  the  value  of  the  flap 
deflection  it  was  computed  for.  By  doing  this  for  flap  de¬ 
flections  from  0"  to  13®  at  an  angle  of  attack  of  0°,  the 
Cjjjac  versus  n  figures  in  Appendix  H  were  obtained.  The 
experimental  points  for  a  particular  E  was  then  found 

by  measuring  the  slope  of  the  corresponding  Cj^gc  versus  n 
figure. 
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Appendix  H 

Experimental  Data  Plots 

This  Appendix  contains  the  experimental  plots  of  the  lift 
coefficient  versus  angle  of  attack  and  the  moment  coefficient 
about  the  aerodynamic  center  versus  flap  deflection  for  all 
seven  models.  The  plots  were  used  to  determine  the  change  in 
zero- lift  angle  of  attack  divided  by  flap  deflection  versus 
flap-chord  ratio  and  the  change  in  the  moment  coefficient 
about  the  aerodynamic  center  divided  by  flap  deflection  versus 
flap-chord  ratio.  The  points  for  the  plots  came  from  reducing 
the  experimental  data  in  Appendix  F  with  the  method  of  reduct¬ 
ion  described  in  Appendix  G. 
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Fig.  22b  Experimental  C^nac  versus  n  plots  used  to  find 
ACmac/n  for  E  -=  ,2 


83 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS-1963-A 


Fig.  2da  Experimental  Cj^  versus  a  plots  used  to  find 
Aou^q/^  for  E  =*  .3 


Fig.  24a  Experimental  C«  versus  a  plots  used  to  find 

for  E  «  .4 


,  Kio  °  -2-75°  for  0”;  td  4j*Ti  ' 
Aal ^  !»  ‘•3. '05®  fori  4*.  to  8®ti| 
Ao,ni-i-5.8*;  for  i0»ito  8rn  i 


Fig.  25a  Experimental  Cj^  versus  a  plots  used  to  find 
Ao^/n  for  E  ■  .5 
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Fig.  27b  Experimental  Cgmc  versus  n  plots  used  to  find 
A<W  /n  for  E  -  .7 


Fig.  28a  Experimental  versus  a  plots  used  to  find 

E  -  .8 
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Fig.  28b  Experimental  versus  ti  plots  used  to  find 

ACjnac/’i  E  -  .8 
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foil  was  chosen  for  this  study.  Results  of  this  study  show 
that  thin-airfoil  theory  does  predict  the  trend  of  the  data 
The  Aolo/'I  versus  E  prediction  favorably  matched  the  data. 
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ed  results.  This  difference  might  be  explained  in  part  by 
the  experimental  setup  used  in  this  study. 
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